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ABSTRACT 


Accurately forecasting future personnel inventory levels by rank and occupational 
specialty is a fundamental prerequisite for development of an effective and functional 
staffing plan. This thesis develops and evaluates univariate time series models to create 
six- and twelve-month forecasts of Marine Corps enlisted manpower levels. Models are 
developed for 44 representative population groups using Holt-Winters exponential 
smoothing, multiplicative decomposition, and Box-Jenkins autoregressive integrated 
moving average (ARIMA) forecasting methods. The forecasts are evaluated against 
actual, out-of-sample inventory levels using several goodness-of-fit indicators including 
Mean Absolute Error Rate (MAPE), Mean Absolute Error (MAE), and Sum of Squared 
Errors (SSE). Among the modeling techniques evaluated, the multiplicative 
decomposition performed the best overall and represents an improvement over the 
Marine Corps’ current naive forecasting method. This thesis recommends Marine Corps 
Systems Command, Total Force Information Technology Systems develop and introduce 
a multiplicative decomposition forecasting model into the Enlisted Staffing Goal Model. 
This forecasting technique should be implemented in phases, starting with the E-1 


through E-4 population groups. 
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EXECUTIVE SUMMARY 


Accurately forecasting future manpower inventory levels in a fundamental 
prerequisite for the development of an effective and functional staffing plan. The 
Enlisted Staffing Goal Model (ESGM) generates a proposed enlisted-force staffing plan 
that mitigates the dilemma of how to optimally meet service-manning requirements 
within the constraints of the Marine Corps’ limited inventory of assignable personnel. In 
order for the staffing plan to be effective, the model’s inputs must be accurate. One of 
these inputs, the forecast inventory available for assignment, is the weakest component of 
the ESGM process. Developing and incorporating a more accurate manpower inventory 
forecasting model will result in a closer to optimal distribution of personnel to meet the 


Marine Corps staffing requirements and improve operational readiness. 


A time series is a collection of observations in time. Although the values of 
individual observations cannot be predicted exactly, the distribution of stochastic time 
series observations commonly follows a discernable pattern. Statistical models can often 
describe these patterns. These models assume that the observations vary randomly about 
an underlying mean value that is a function of time. The time series may also be 
characterized by one or more behavioral components that may be isolated, modeled and 


incorporated into the forecast algorithm. 


Univariate time series forecasting models forecasting models make predictions by 
extrapolating the past behavior of a single variable of interest. Forecasting is appropriate 
for stochastic time series data when the underlying causes of variation do not change 
significantly over time. The forecasting process consists of five steps. This heuristic 
method includes formulating the problem, obtaining data, selecting and applying 
forecasting methods, evaluating models, and using forecasts. Guided by this 


development framework, analysts can produce more accurate and efficient forecasts. 


This thesis develops and evaluates manpower inventory forecasting models for 44 
representative Marine Corps enlisted population groups using three univariate time series 


forecasting techniques. These methods are the Holt-Winters exponential smoothing, 


Xlil 


multiplicative decomposition, and Box-Jenkins autoregressive integrated moving average 
forecasting models. Historical personnel strength levels obtained from the Marine Corps 


Total Force Data Warehouse comprise the observation sets employed in this thesis. 


The Holt-Winters exponential smoothing method is a common univariate time 
series forecasting technique. This model smoothes irregular fluctuations by using 
weighted averages in an exponentially decreasing manner. The technique creates 
forecasts by estimating and extrapolating a linear trend while adjusting the data in each 


period by estimated seasonal indices. 


Multiplicative decomposition is another effective method of forecasting univariate 
time series data. This technique is based on the concept that the underlying factors can 
be identified and isolated. These factors are trend, cycles, seasonality, and random 
variation. The analyst develops a model by first identifying and removing the component 
effects from the data. After these effects are isolated, the analyst creates a forecast by 


reassembling the components. 


The Box-Jenkins technique is a sophisticated approach by which to analyze time 
series data and extrapolate a forecast. This methodology provides a framework for 
preparing data, selecting a model, and creating forecasts commonly referred to as 
autoregressive integrated moving average (ARIMA) forecasting. This technique assumes 
any given time series observation can be modeled as a function of its past values and 


current and past values of random errors. 


Models developed in this study using the Holt-Winters exponential smoothing, 
the multiplicative decomposition, and the Box-Jenkins forecasting techniques are 
compared to the current forecasting method and each other applying a variety of 
statistical goodness-of-fit measurements. The evaluation techniques used are the Mean 
Absolute Percentage Error, the Sum of Squared Errors, and the Mean Absolute Error. 
Among the forecasting techniques evaluated, the multiplicative decomposition method 
performed best overall, and represents an improvement over the Marine Corp’s current 


forecasting method. 
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I. INTRODUCTION 


A. BACKGROUND 
The purpose of the Marine Corps assignment process is to satisfy manpower 
readiness requirements through the best distribution of available personnel. Accurately 
forecasting future personnel inventory levels by rank and occupational specialty is a 
fundamental prerequisite for development of an effective and functional staffing plan. 
Manpower forecasting seeks to increase an organization’s human resource 


options and reduce the penalty costs of lost opportunity and lowered 
performance in managing those resources (Burack & Mathys, 1980). 


The Enlisted Staffing Goal Model (ESGM) generates a proposed enlisted-force 
staffing plan that mitigates the dilemma of how to optimally meet service-manning 
requirements within the constraints of the Marine Corps’ limited assignable personnel 
inventory. In order for the staffing plan to be effective, the model’s inputs must be 
accurate. These inputs include manpower requirements, staffing policies, and projected 


available personnel inventory. 


The ESGM’s forecasted manpower inventory is the model’s leading weakness. 
Currently this supply is represented by a snapshot of the existing personnel inventory 
captured six-months prior to the execution phase of the assignment process. This naive 
forecasting method does not consider the myriad of factors that affect personnel 
inventory levels over time including promotions, reductions, reenlistments, extensions, 
deaths, MOS changes, and administrative discharges. This negatively impacts the utility 
of the staffing plan. 


Without an accurate inventory-forecasting model, the ESGM currently produces 
unrealistic staffing goals that are frequently unachievable at the time of execution. This 
shortcoming significantly impairs the realization of the model’s primary objective — 
ameliorating the manpower system’s mismatch between requirements and the assignable 
population. Developing and incorporating a more accurate inventory forecasting method 
will result in a closer to optimal distribution of personnel to meet the Marine Corps’ 


manpower requirements and improve operational readiness. 


B. PURPOSE 

The purpose of this thesis is to improve the output of the Marine Corps Enlisted 
Staffing Goal Model by more accurately predicting the personnel resources available for 
assignment to future manpower requirements. The primary research question is the 
following: What univariate time series modeling technique is the most effective method 
to forecast Marine Corps enlisted personnel inventory levels by rank and occupational 
specialty, six and twelve months in the future? 
C. SCOPE AND METHODOLGY 

This thesis develops and evaluates manpower inventory forecasting models for 44 
representative Marine Corps enlisted population groups using three univariate time series 
forecasting models. These techniques are the Holt-Winters exponential smoothing, 
multiplicative decomposition, and Box-Jenkins autoregressive integrated moving average 
forecasting methods. Historical personnel strength levels obtained from the Marine 
Corps Total Force Data Warehouse comprise the observation sets employed in this study. 
These time series contain 48 monthly observations starting in October 2001 and 
concluding in September 2005. The forecasts are evaluated against actual, out-of-sample 
inventory levels using several goodness-of-fit measures including Mean Absolute Error 
Rate, Sum of Squared Errors, and Mean Absolute Error. 
D. BENEFITS OF THE STUDY 

By developing and evaluating various manpower forecasting models, this study 
will provide the necessary information required to integrate a more effective enlisted 
personnel strength prediction into the Enlisted Staffing Goal Model. Incorporating this 
forecasting model will result in a closer to optimal distribution of the assignable 
inventory to meet the Marine Corps’ manpower requirements and improve readiness. 
E. ORGANIZATION 

This introductory chapter provides the reader with the motivation, objective, and 
organization of this thesis. Chapter II offers a synopsis of time series analysis and the 
forecast model development process. Chapter III presents an overview of the Marine 
Corps Enlisted Staffing Goal Model. Chapter IV introduces the data used to develop the 
study’s forecasting models. Chapters V through VII discuss the Holt-Winters 


exponential smoothing, multiplicative decomposition, and Box-Jenkins modeling 
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techniques, respectively, and demonstrate forecast development. Chapter VHI reviews 
model evaluation techniques and presents a comparison and analysis of the forecast 
results. The thesis concludes with recommendations and potential areas for further 


research. 
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I. TIME-SERIES ANALYSIS AND FORECASTING 


A. OVERVIEW 

A time series is a sequence of observations obtained through measurements often 
recorded at equally spaced intervals. Often, time series data have characteristics that 
facilitate forecasting. These include seasonality, underlying trends, and relationships 
with past observations or other causal variables. Analysts can improve time series 
forecasts if they understand the nature of these components and identify the model that 


will best exploit the data’s characteristics. 


The purpose of this chapter is to provide a synopsis of time series analysis and 
forecasting. The first section discusses the characteristics of time series data. It reviews 
the common components useful in creating effective forecasts such as trend, seasonality, 
cyclical behavior, and irregular fluctuations. The chapter concludes with an introduction 
to time series forecasting and an overview of the forecasting model development process. 
B. CHARACTERISTICS OF TIME-SERIES DATA 

A time series is a “collection of observations made sequentially in time” 
(Chatfield, 1996). Examples are records of local daily rainfall levels, the quarterly U.S. 
Gross Domestic Product, and the monthly Marine Corps personnel strength for a 
particular rank and MOS. Time series analysis provides tools for choosing a 


representative model and producing forecasts. 


There are two kinds of time series data: 


@ Continuous, where the data contain an observation at every instant of time, 
e.g., seismic activity recorded on a seismogram. 


° Discrete, where the data contain observations taken at intervals, e.g., 
monthly crime figures. 


Unless the data are purely random, observations in a time series are normally 
correlated and successive observations may be partly determined by past values 
(Chatfield, 1996). For example, the meteorological factors that affect the temperature on 
any given day are likely to exert some influence on the following day’s weather. Thus, 


historical temperature observations are beneficial in forecasting future temperatures. 


A time series is deterministic if it contains no random or probabilistic 
characteristics but proceeds in a fixed, predictable fashion (Chatfield, 1996). An example 
of a deterministic time series would be the data collected while conducting a classical 
physics experiment such as one demonstrating Newton’s law of motion (Gujarati, 2003). 
More applicable to econometric applications are stochastic time series. Stochastic 
variables have indeterminate or random aspects. Although the values of individual 
observations cannot be predicted exactly, measuring the distribution of the observations 
may follow a predictable pattern. Statistical models can describe these patterns. These 
models assume that observations vary randomly about an underlying mean value that is a 
function of time. Time series data can also be characterized by one or more behavioral 
components: trend, seasonality, cyclical behavior, and random noise. 

1. Trend Component 

Trend is the general drift or tendency observed in a set of data over time. It is the 
underlying direction (an upward or downward tendency) and the degree of change in an 
observation set when consideration has been made for other components. Graphing a 
time series can be a useful and simple method of identifying the trend of a particular data 
set. Figure | indicates an upward trend of the U.S. Gross Domestic Product over a ten- 
year span. Analysts can also discern trends by dividing the data set into a number of 
ranges, and calculating the mean for each span. A consistent increase or decrease in the 


mean for the successive ranges indicates trend. 
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Figure 1. | Annual U.S. Gross Domestic Product ($ billions) 1984-2004 
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Trends in business or economic series may be due to a growth or contraction 
process. Trends in service manpower levels may be attributed to external economic 
factors or shifts in policy due to technical innovation, downsizing, or an increased or 
decreased operational requirement for certain occupational specialties. 

2. Seasonal Component 

In time series data, the seasonal component is the element of variation in a data 
set that is dependent on the time of year. Seasonality is quite common in econometric 
time series. It is less common in engineering and scientific data. This component recurs 
annually, with possible variations in amplitude. Seasonality is attributable to the change 
of seasons and/or the timing of such events as holidays or the start or completion of the 
school term. For example, the cost of fresh produce, retail sales levels, average daily 
rainfall amounts, and unemployment figures all demonstrate seasonal variation. Figure 2 
effectively illustrates the effect of seasonality on the California unemployment rate over a 


twenty-year period. 
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Figure 2. Monthly California Unemployment Rate 


Incorporating seasonality in a forecast is useful when the time series has a 
discernible seasonal component. When the data contain a seasonal effect, it is useful to 
separate the seasonality from the other components in the time series. This enables the 
analyst to estimate and account for seasonal patterns. 

BE Cyclical Component 

Cyclical behavior describes any non-seasonal component that oscillates in a 
recognizable pattern. The 11-year sunspot cycle has been long recognized as naturally 
occurring cyclical activity. More ambiguous is the 5 to 7 year business cycle that a 
number of economists hypothesize influence global economic activity. If the data 
include a discernible cyclical component, the time series should span enough cycles to 
accurately model and forecast its effects (Yaffee, 2000). Cyclical behavior in which the 
oscillations extend over a very long period (such as 20 years) can often be accurately 
modeled as a trend for short-term forecasts (Chatfield, 1996). 

4. Irregular Component 

A significant component of stochastic time series data is irregular fluctuation. 
This random noise is what remains after the other components of the series (trend, 
seasonality, and cyclical behavior) are estimated and eliminated. It results from 
fluctuations in the series that are neither systematic nor predictable. While the irregular 
component typically has a modest impact on attempts to analyze, model and forecast time 
series data, it can sometimes have a significant effect. The effect of the 1973 OPEC 
(Organization of Petroleum Exporting Countries) oil embargo on the U.S. economy is an 
example of the substantial consequences irregular fluctuations can occasionally inflict on 
time series data. 

C. TIME-SERIES FORECASTING 

Univariate time series forecasting models make predictions by extrapolating the 
past behavior of the values of a particular single variable of interest (Moore & 
Weatherford, 2001). Successive observations in econometric time series are normally 
not independent and predictions may be made from previous observations (Chatfield, 
1996). While exact forecasts are possible with deterministic time series, forecasts of 
stochastic time series are limited to “conditional statements about the future based on 


specific assumptions” (Chatfield, 1996). According to Armstrong (2001), “the basic 


assumption is that the variable will continue in the future as it has behaved in the past.” 
Specifically, time series forecasts are appropriate for stochastic data where the underlying 
causes of variation — trend, cyclical behavior, seasonality, and irregular fluctuations — do 
not change significantly in time (Jenson, 2004). Hence, modeling is often more 
appropriate for short-term than for long-term forecasting. 

1. The Forecasting Model Development Process 

In his book Forecasting Principles, Armstrong (2001) proposes a_ heuristic 
forecasting process. Guided by the five general steps of Armstrong’s development 
framework, analysts can produce more accurate and efficient forecasts. These steps are 
graphically depicted in Figure 3 and include formulating the problem, obtaining 
information, selecting and applying forecasting methods, evaluating models, and using 
the forecasts. The remainder of this chapter summarizes Armstrong’s forecasting 


principles as they apply to this thesis. 
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Figure 3. | Forecast Development Process (From: Armstrong, 2001) 


2. Formulate the Problem 

The first step in Armstrong’s (2001) forecast development process is to formulate 
the problem. The analyst must clarify the user’s objectives for the forecast and the 
decisions associated with these particular goals. Of particular note is the initial 
assessment as to whether a practical forecast is possible for a particular problem. 
Forecasting should be avoided if previous studies indicate that the scenario is unlikely to 


yield a beneficial result. Additionally the analyst should structure the problem 
, 


appropriately in order to forecast effectively and produce results in a format constructive 
for the user. This might include decomposing the underlying factors affecting the 
situation, identification of significant relationships between causal variables, and 
recognition of trends or other components of time series data. 

3. Obtaining Information 

After appropriately formulating the problem, the next step in Armstrong’s (2001) 
forecast development process is the identification, collection, and preparation of data. 
The analyst can employ theory and previous research to identify potential explanatory 
variables. These data should correspond to the forecasting problem and be objective in 
nature. Often the best predictor of future behavior, given a particular scenario, is past 
behavior. Once a source is identified, the forecaster must collect all pertinent, suitable 
and reliable data and prepare it for the forecasting process. Preparing the data for 
modeling includes scrubbing the series, adjusting for errors and missing observations 
where required and employing statistical techniques to adjust for unsystematic past 
events. A hurricane would be an example of an unsystematic past event when analyzing 
regional retail sales time series data. Some models may also require data transformation 
in order to render the series stationary in mean, variance, and autocovariance (Yaffee, 
2000). 

4. Select and Implement Forecasting Methods 

Often many methods are applicable to a particular forecasting situation. The 
forecaster should list and consider all of the selection criteria before selecting a 
forecasting method. Although accuracy is an essential measure of effectiveness, other 
criteria such as cost, development time, expertise available, ease of use, and ease of 
implementation are important considerations as well. Quantitative forecasting methods 
tend to be less subjective than qualitative techniques and are suitable for occasions when 
germane data are readily available. If feasible, the analyst should evaluate the forecasts 
of several methods under comparable conditions noting the accuracy, cost, and usability 


of each technique. 


Occam’s Razor states that given two equally predictive theories, an individual 


should choose the simplest. Thus, to improve accuracy, assist user comprehension, 
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decrease errors, and reduce expenses, the forecasting model should be as uncomplicated 


as practical — using few variables and simple relationships between among them. 


Additionally, the model should be tailored to the requisite planning horizon. 
Short-term forecasts should place added weight on the most recent observations. Distant 
forecast lengths should emphasize applicable long-term trends and cycles. Although 
there are no established guidelines regarding minimum sample size, the series should 
contain sufficient observations for accurate parameter estimation. If the series contains 
seasonality or other cycles, the data should span enough oscillations to effectively model 
the cyclical behavior (Yaffee, 2000). Thus seasonal or cyclical processes “require more 
observations than non-seasonal ones” (Yaffee, 2000). In addition to Armstrong’s (2001) 
criteria for model selection and implementation above, Chatfield (1996) recommends 
choosing methods that account for the number of series to be forecast and the resources 
required per forecast, suggesting that “univariate forecasts are particularly suitable when 
there are large numbers of series to be forecast.” 

3s Evaluate Methods 

Methods for evaluating forecast models are based on conventional analytical 
procedures. Armstrong (2001) suggests evaluating at least two “reasonable” forecasting 
techniques including the current method if applicable. The analyst should evaluate the 
models under situations that match the forecasting problem. One effective method is to 
back-forecast employing older observations in the time series comparing the forecast fit 
against the actual values of recent, out-of-sample observations (Wooldridge, 2003). 
When evaluating the effectiveness of forecast methods across several series, the analyst 
should utilize error measurements that consider the scale of the observations (Armstrong, 
2001). This measurement must not be subject to distortion due to a disparity in series 
observation value magnitudes (thus weighting some samples more heavily). <A 


measurement that incorporates percentage error is appropriate for these circumstances. 


If the forecasting methods evaluated are not adequate for the problem, the analyst 
should consider another technique. Finally, the analyst and the end user should consider 
the costs and benefits of each model evaluated in order to select the most appropriate 


method. 
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6. Use Forecasts 

After selecting a model, the analyst should present the forecast and the supporting 
data in a straightforward and comprehensible format, such as a graph. The 
documentation should provide a complete and clear explanation of the data, methods, and 
assumptions. The forecaster should review the forecast accuracy periodically, updating 
the parameters of the model as the data or situation changes and actively seek user 


feedback. 
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Hl. THE ENLISTED STAFFING GOAL MODEL 


A. OVERVIEW 

The Enlisted Staffing Goal Model (ESGM) generates a proposed enlisted force- 
staffing plan that resolves the dilemma of how to optimally meet service-manning 
requirements within the constraints of the Marine Corps’ limited assignable personnel 
inventory. The primary focus of this chapter is to provide an overview of the ESGM. 
First, the model’s purpose and utility are examined. Next, the ESGM is examined from a 
system perspective, including input, processes and output. Finally, the chapter concludes 
with an examination of effectiveness of the model and identification of a significant 
limitation that impedes the ESGM’s value. 
B. FUNCTIONALITY AND USE 

The purpose of the Marine Corps’ inventory assignment process is to make the 
best distribution of current assignable inventory to meet the Marine Corps authorized 
strength requirement according to current staffing precedence policies. The foundations 
of personnel requirements are established in the individual unit Tables of Organization 
(T/O’s). These documents are based upon each unit’s defined mission. The sum of the 
service’s T/O’s typically exceed the number of personnel authorized by Congress and 
available for assignment. Currently, the Marine Corps has over 177,000 active duty 
service members. Of these, approximately 144,000 Marines are available to staff 
154,000 structured billet requirements. The manning plan developed annually by 
Headquarters Marine Corps is called the Authorized Strength Report (ASR) and is 
developed based on T/O’s, Congressional authorization, operational requirements and 
service policies. The quantities and skills of the existing inventory of Marines are not 
considered when producing the ASR. Thus there is routinely a mismatch between the 
actual number of billets authorized and the assignable population. The Manpower and 
Systems Support Section of the Headquarters Marine Corps’ Enlisted Assignments 
Branch employs the ESGM to manage the predicament of how to allocate the assignable 
personnel inventory to the billet requirements with the best skill fit (as defined by rank 


and military occupational specialty (MOS)). 


13 


The purpose of the ESGM is to produce a set of targets called staffing goals, 
which represent a billet requirement that can be filled using the available inventory. This 
staffing plan is provided to subordinate organizations as guidance from Headquarters 
Marine Corps indicating which billets the Enlisted Assignments Branch intends to fill at 
each unit annually. Staffing goals are perceived as a commitment and frequently 
described metaphorically as a “debt” which the Enlisted Assignment branch “repays” to 
the operating force and supporting commands each staffing period using the assignable 
inventory as “currency” (HQMC, MMEA-S, 2005). 

C. ESGM STRUCTURE AND PROCESSING OVERVIEW 

The staffing goals also guide the individual detailers in the execution of the 
assignment. However, in order for the staffing goals to be practical and functional, the 
inputs must be effective and accurate. The ESGM’s inputs include requirements, 


assignment and staffing policies, and projected available personnel inventory. 


Manpower requirements input is provided in the form of the ASR which 
constrains the combat-based, Marine Corps T/O personnel requirements by (a) the 
manpower levels authorized by Congress and (b) the “Transients, Trainees, Prisoners and 
Patients” (T2P2) account. T2P2 estimates the number of Marines unavailable for 
assignment annually for time spent executing orders, training in excess of 20 weeks, 
hospitalization for greater than 30 days, or incarceration for between 30 and 180 days. 


The ASR identifies billet requirements by unit, pay grade, MOS, and number authorized. 


The model also integrates assignment and staffing policies to balance the 
insufficient supply with the demand. These are a set of business rules and constraints that 
optimize the “fill” and “fit” in accordance with the Commandant’s priorities and billet 
requirements. The “fill” is determined by the staffing precedence order (Marine Corps 
Order 5320.12E, Precedence Levels for Manning and Staffing) that specifies which 
particular units will be staffed at 100% of T/O level (Excepted Commands), which units 
will be staffed at 95% of T/O level (Priority Commands), and which units will be staffed 
with the remaining inventory on a proportionate share basis (Pro-Share Commands). The 
“fit” is determined by a set of restrictions and substitution rules, which may permit grade 


and MOS substitutions or establish gender restrictions for particular billets. 
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The third input, the forecasted assignable manpower inventory, is provided in the 
form of an extract of personnel information from the Marine Corps Total Force System 
Operational Data Store Enterprise (MCTFS/ODSE). The ESGM processes the data and 
creates a forecasted inventory of assignable personnel based on the specific Marines’ 
remaining service obligations. Specifically, the model considers an individual eligible to 
fill a staffing goal for his or her current rank and MOS if his or her end of active service 


(EAS) date exceeds the current date by six months or more. 


The optimization process has two objectives: maximize the fill of billets and 
maximize the fit of individuals to the billets. The ESGM maximizes the fill while 
honoring the policy-defined priority. Billets with higher priority are filled before lower 
priority billets. Within any priority level, shortages are shared in proportion to the target 
if allowed by “fit” rules. The result of this process is a set of feasible quotas, which the 
ESGM then uses to maximize the “fit” of individuals. The ESGM endeavors to allocate 
individuals within the fill constraints to the highest level of desirability as measured by 
occupational specialty. Finally, the model attempts to allocate individuals with a perfect 


grade match proportionately to the targets. 


Included among the outputs generated by the ESGM are a database containing all 
information required to run the Enlisted Assignment Model (EAM), data posted to the 
Monitor Assignment Support System (MASS), and the Command Staffing Report (CSR). 
The CSR is a critical document used by both the assignment monitors and individual 
commands. It specifies which T/O billets are authorized and which billets received a 
staffing goal for a particular staffing period. 

D. ACCURACY 

The ESGM’s forecasted manpower inventory is a major weakness of the model. 
Currently this supply is represented by a snapshot of the existing personnel inventory six 
months prior to the execution phase of the assignment process. Other than removing 
those Marines with an approved separation date, the model does not forecast the 
projected inventory at the time of execution. This naive forecasting method does not 
consider the myriad of factors that affect personnel inventory levels over time including 
promotions, reductions, reenlistments, extensions, deaths, MOS changes, and 


administrative discharges. This affects the utility of the staffing plan. 
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Without an effectual inventory-forecasting model, the ESGM currently produces 
unrealistic staffing goals that are unachievable at the time of execution. This 
shortcoming significantly impairs the realization of the model’s primary objective — 
ameliorating the manpower system’s mismatch between requirements and the assignable 
population. Developing and incorporating a more accurate inventory forecasting method 
will result in a closer to optimal distribution of personnel to meet the Marine Corps’ 


manpower requirements. 


16 


IV. DESCRIPTION OF DATA 


A. MONTHLY MANPOWER INVENTORY LEVELS 

This thesis develops and evaluates models for the 44 MOS/paygrade 
combinations detailed in Table 1. The Marine Corps Enlisted Assignments Branch 
considers these groups a representative sample of the overall enlisted Marine Corps 
population. The Appendix summarizes the time series forecasts constructed for each 
population. Monthly historical personnel inventory level time series were obtained from 


the Marine Corps Total Force Data Warehouse (TFDW). 


Marine Corps Manpower 
Inventory Time Series Analyzed 


E5 0211 ¢ E7 0369 ¢ E8 6019 

E6 0211 ¢ E8 0369 ¢ E9 6019 

E7 0211 ¢ E9 0369 ¢ E1-E3 6092 
E8 0211 ¢ =E5 2336 ¢ E4 6092 
E1-E3 0311 ¢ E6 2336 ¢ E5 6092 

E4 0311 ¢ E7 2336 ¢ E6 6092 

E5 0311 ¢ E8 2336 ¢ E7 6092 
E1-E3 0321 ¢ E9 2336 ¢ E1-E3 7257 
E4 0321 ¢ E1-EK3 5811 ¢ E4 7257 

E5 0321 ¢ E45811 e = E5 7257 

E6 0321 e E55811 ¢ E6 7257 

E7 0321 ¢ EK65811 ¢ E7 7257 

E8 0321 ¢ E75811 ¢ E8 7291 

E9 0321 ¢ E8 5811 ¢ E9 7291 

E6 0369 ¢ E95811 

Table 1. MOS/Paygrade Population Inventory Time Series Analyzed 


Each time series contains monthly manpower inventory levels with 48 
observations starting in October 2001 and concluding in September 2005. The modeling 
techniques evaluated in this thesis employ the first 36 observations to develop forecasts. 
The final 12 observations in each series are reserved for evaluating each model’s forecast 
accuracy with respect to actual inventory levels. The period includes a 6-month interval 


in which the Marine Corps enacted a stop-loss policy in response to combat operations in 
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Iraq. This policy involuntarily extended Marines on active duty, substantially increasing 
the inventory levels for several population groups from February 2003 until July 2003. 
B. E-5 0311 TIME SERIES 

The following chapters refer to the 0311 E-5 (Rifleman sergeant) time series in 
order to demonstrate the development of the exponential smoothing, multiplicative 
decomposition, and Box-Jenkins forecasting models. Table 2 provides the monthly 
inventory levels for the 0311 E-5 population group retrieved from the TFDW. ‘The 
graphical representation of the data provided in Figure 5 demonstrates the volatility of the 
inventory levels and highlights the time series’ significant irregular component resulting 


from stop-loss. 


[Period __|Month __[Inventory|Period [Month _|Inventory| 
| 63/1/02] 1857] 30] 3/1/04] 1889) 
| 85/1/02] 1877] 82] 5/1/04] 1898) 
| 6/1/02] 1866] 33] 6/1/04] 1856) 


Table 2. USMC 0311 E-5 Monthly End Strength, Oct 01 — Sept 05 
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Figure 4. 0311 E-5 Inventory Levels Oct 01 — Sept 05 
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V. EXPONENTIAL SMOOTHING FORECASTING 


A. OVERVIEW 

Exponential smoothing is a common univariate time series forecasting technique. 
This method smoothes irregular fluctuations by using weighted averages of time series 
observations in an exponentially decreasing manner. The older an observation is, the less 
influence it has on the forecast. More advanced exponential smoothing models endeavor 
to identify and model time series components and incorporate them into the forecast 
(Yaffee, 2001). Exponential smoothing is appropriate for automatic, short-term 
forecasting of frequently used data where the underlying causes of variation do not 


change markedly over time. 


This chapter reviews exponential smoothing forecasting techniques. An 
examination of the simple exponential smoothing model provides the technical 
foundation for the remainder of the chapter. The following sections discuss linear and 
seasonal exponential smoothing methods. The chapter concludes by demonstrating the 
application of the exponential smoothing forecasting process. It employs the Holt- 
Winters technique to forecast manpower inventory levels for the E-5, 0311 Marine Corps 
enlisted population group introduced in Chapter IV. 

B. SIMPLE EXPONENTIAL SMOOTHING 

Exponential smoothing in its most basic form is appropriate for forecasting time 
series with no seasonality or trend components (Chatfield, 1996). It provides a starting 
point for the more advanced smoothing techniques evaluated in this thesis. Exponential 
smoothing places more weight on the most recent observations and exponentially less 
weight on older values. To examine this method in detail, consider a time series 


{X,,X,,...,x,} where x, is the current observation. The analyst can extrapolate the value 


at one time period in the future, x,,,, 


as a weighted sum of the previous observations. 


Consider the formula: 





Kg = COX FGM FO Ky Fev. (5.1) 
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where c, are weights. In Equation (5.1), the weights decay exponentially from the most 


recent to the most distant observation. These weights sum to one, and the proportion of 
the latest observation taken is the constant @, commonly referred to as the damping 


factor or the smoothing constant: 


c, =a(1—-a) 
where: i=1,2,..., (5.2) 
and: O<a<l. 


Thus, Equation (5.1) becomes: 


Sip = aX, ¢a(l-a)x,,+a(l-a@)’x,, 4+... (5.3) 


This formula presumes an unlimited number of observations when in actuality the 
forecast requires only a finite number (Chatfield, 1996). Therefore, Equation (5.3) is 


restated in the recurrence form as: 


K,4, = ax, +(l-@)lax,_,+aU—-a)x,_,+...] (5.4) 
=ax,+(1-@)x,_. , 


This equation requires an estimate of x,. A straightforward and effective option 


is X,=x,. Equation (5.4) can be rewritten to incorporate error correction: 
X= %, + a(x, —%,) (5.5) 


where, beginning with t =1, the forecast for the next period, * 


t+1? 


equals the sum of the 
forecast for the current time period, x,, and the weighted error of the current period’s 
forecast, a@(x,—X,). Forecasts are computed in this manner for all observations in the 


time series. 


The choice of the smoothing constant, @ , impacts the accuracy of the forecast 
significantly. Alpha is calculated to minimize the sum of discrepancies between the time 
series’s actual observations and the forecast values. This computation takes the form of 


an optimization program: 


Ze 


t 
Minimize : ye (x, -£,)? 


i=l 


Subject to: X,,, =X, +a(x,-%,), (5.6) 
i=1,2,...,¢-1, 
and O<sasl. 


C, HOLT’S LINEAR EXPONENTIAL SMOOTHING 

Simple exponential smoothing does not account for the trend component common 
in many time series. By estimating and extrapolating the series’ linear trend, the Holt 
linear exponential smoothing method often produces a forecast more accurate than simple 
exponential smoothing. The model’s intercept and slope are calculated from weighted 
averages of previous observations, using separate smoothing constants for each. These 


constants are denoted by @ and f respectively. The equations required to calculate the 


level, L,, and linear trend, T,, are: 


LL, =ax,+(-ay(L,+T,) (5.7) 

where: 0<a<l 

T= f(44-L)+d-AFT, 

t+1 B( t+1 ) ( B) t (5.8) 

where: 0< 6 <1. 

Initial estimates are required for L, and 7,. Simple options are: 

Lae (5.9) 
= 5 hs (5.10) 


Thus, the forecast for the next time period is calculated using Equation (5.11), 


below. Predictions for i time periods in the future are forecast with Equation (5.12). 


kK =L,,4+T 


t+1 ~~ 7441 +1 


(5.11) 


K,4 =£,,+ OF, (5.12) 


+1 


As with the simple exponential smoothing model, the values of the damping 
factors, a and f, are important determinates of the model’s accuracy. Minimizing the 
disparities between the model’s actual and predicted values will optimize these constants. 


Equation (5.13) summarizes this optimization program. 
Pa 


Minimize: er -%) 
i=] 

Subject to: L,,, =ax,+(—a)(L,+T), 

T= B(L,,-L)+d-AE, 

Ki =L,,+Ti5 

i=1,2,...,t-l, 

O0<a<l, 


and 0<f<1 
D. HOLT-WINTERS SEASONAL EXPONENTIAL SMOOTHING 


The Holt-Winters forecasting technique is an extension of linear exponential 


(5.13) 


smoothing which accommodates a seasonal component present in a time series. This 
model forecasts by estimating and extrapolating a linear trend while adjusting the data in 
each period by estimated seasonal indices. The intercept, slope and seasonality are 
computed using weighted averages of previous observation values, applying separate 


smoothing constants to each. These parameters are denoted a, (, and y respectively. 


Level, trend and seasonality are designated L,, T, 


;, and S, at time ¢, while the 

subscript s is the periodicity of the seasonality. Therefore, if the time series consists of 

monthly observations, s = 12. If the series is comprised of quarterly observations, s = 4. 

Equation (5.7) is modified to accommodate seasonality in the following manner!: 
LL, =a(x,/S,,)+d-a)(L,+T,) 


where: 0<a<\l. 


(5.14) 


This calculation requires the extrapolation of T,, L, and S, for the first complete 


cycle of seasonality — the first twelve observations for monthly time series data. 


Equations (5.15) (5.16) and (5.17) provide straightforward methods of estimation. 


f= 457%; (5.15) 

jamb tes, (5.16) 
Ss 

Pig = X53 Lies (5.17) 


| This thesis describes the multiplicative forecasting method. For a description of the additive method, 
see Yaffee (2001). 
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where :i =1,2,...,5 
Trend is computed in the same fashion as in linear exponential smoothing: 


T= BE —h)+ 0 AT, 


(5.18) 
where: 0< #<1. 


After extrapolating the values, the analyst can compute the forecast for the next 


time period: 


ae = (L, + T, Sra, 7 (5. 19) 
Seasonality is updated after observing the actual value of x;41: 


Sra a V(X ‘ E41) a Ura vere 


(5.20) 
where:0<y<l. 


Forecasts for time periods beyond the next are calculated in the following manner: 


£,, =(L,,, + 2T.,,)S 


t+1 t+ t+2-s 
443 = (L., a OTe) Oe 
Kies = (L,., F sT,,,)S, (5.21) 
Nica = (Ga a (s ae DT... Svat 


Asics raz (L,., 7 (s as 2) )Si49-s 


The smoothing constants are calculated via an optimization function similar to the 


simple and linear exponential smoothing models: 


Minimize : ¥e" -%) 
i=l 

Subject to: L,,, =a(x,/S, .)+(U-ay(L, +T,), 
T= BCL, -L)+0d-AE, 
Sin = 1% LS 
Ki = La +b )Spt-5> (5.22) 
i=st+lLs+4+2,..t-l, 
O0<a<l, 
0<f<l, 

and O<y<l. 


t+l-s? 


pes, 


E. MANPOWER FORECASTS WITH HOLT-WINTERS 

Holt-Winters exponential smoothing techniques are applicable to forecasting 
manpower inventory levels for many Marine Corps enlisted population groups. This 
section demonstrates the process to predict E-5 0311 inventory levels six and twelve 


months in the future. 


This time series contains monthly observations; therefore the periodicity of 
seasonality is twelve, i.e. s=12. Only the first 36 observations will be used to create the 
forecast. The final twelve observations will be reserved to evaluate the forecast accuracy 
and are considered “out of sample” for the model. Therefore, for the purposes of the 
forecast, the most recent observation is period 36, the inventory level on September 1, 


2004. Consequently, t= 36 and the current observation is denoted x,,. 


Level, trend, and seasonality for the first twelve observations, {x,,x,,...,%,,}, are 


computed using Equations (5.9), (5.16), and (5.17). These figures are presented in Table 
3. 


e857, 1857] 45| 1.0000] 


-__e[_i877|__1877|__-11] 1.0000 
__3|_1866[ 1866] _-57| 1.0000 





Table 3. Level, Trend & Seasonality Calculations; {x,,x,,...,%,,} 
1,13 and 7T)3 are computed after observing x;2 using Equations (5.14) and (5.18). 
At this point, the analyst has sufficient information to estimate X,,, which is the forecast 
inventory level for October 1, 2002. Using Equation (5.19), the estimate, x,,, is given by 
(113+7}3)*S1. After observing the actual value of x13, the forecaster must update the 


seasonal effect for October, $13, with Equation (5.20). This process is repeated for each 
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remaining observation in the time series, {X,;,%,4,.-.,%,,}. These calculations are 


presented in Table 4. 





Table 4. Level, Trend, Seasonality & Forecast Calculations; {x,,,%,4,.--. X36} 


The forecaster should now optimize the smoothing constants by minimizing the 
disparities between the model’s predicted values and the time series’ actual observations. 
The optimization function specified in Equation (5.22) is easily accomplished with the 
Microsoft Excel Solver plug-in. The forecasting problem’s optimized smoothing 
constants are presented in Table 5. The small value for gamma is indicative of the time 


series’s nominal seasonal component. 


Alpha | 0.9774 
0.46994 





IGamma | 0.08733 


Table 5. Optimized Smoothing Constants 


Using Equation (5.21), the analyst can now produce forecasts. The formulas for 


six and twelve month forecasts are presented as Equations (5.23) and (5.24) respectively, 
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below. Table 6 presents the predicted inventory level for each forecast period, the 
observed (actual) inventory level, and the absolute percentage error. The absolute 
percentage error, described in Equation (5.25), is one of many possible measures of 
forecast effectiveness. Figure 6 provides a graphical presentation of the observed and 


forecast inventory levels for each period. 


¥qq = (Lay + 6T 57) S39 (5.23) 
Kyg = (UL +127, ) Sy, (5.24) 
\(~,- 4) 
APE =———— * 100 (3.25) 
Xx, 


L 


Holt Winters Exponential Smoothing Forecast Results 
Forecast Period |Actual Inventory |Predicted Inventory|Absolute Percentage Error 


6-Month 1869 1762 5.72% 
12-Month 1880 1667 11.33% 


Table 6. E-5 0311 Holt Winter’s Exponential Smoothing Forecast Results 





F. CHAPTER SUMMARY 

Exponential smoothing is a common technique used to extrapolate a forecast from 
a historical time series. This method considers the entire past in the model, but weighs 
recent observations more heavily than older observations. Simple exponential smoothing 
models are appropriate for time series with nominal trend and seasonality components. 
Holt’s linear smoothing technique isolates and models the series’s trend, often producing 
a more accurate forecast than simple exponential smoothing. The Holt-Winters method 
is an extension of linear exponential smoothing which accommodates a seasonal 
component present in a time series. These models are appropriate for automatic, short- 
term forecasting of frequently used data where the underlying causes of variation are not 


changing markedly in time. 
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VI. TIME SERIES DECOMPOSITION FORECASTING 


A. OVERVIEW 

Time series decomposition is a forecasting and analysis method frequently 
employed by the U.S. Bureaus of Census and Labor Statistics. Frederick R. Macaulay 
pioneered this technique in the 1920’s at the National Bureau of Economic Research 
(Yaffee, 2000). The purpose of the research was the seasonal adjustment of time series 
data for analysis; however, refinements to the work have yielded practical and effective 


techniques to analyze cyclical behavior and create univariate forecasts. 


The fundamental assumption of the decomposition technique is that every time 
series consists of a number of component parts that are related in some manner. As 
introduced in Chapter II, these components are trend (7), cycle (C), seasonality (S$), and 
irregular (J) or random variation. The decomposition methodology views a time series as 


the sum or product of these attributes:2 
x=P*S*C#T (6.1) 


This chapter examines time series decomposition forecasting techniques. Using 
the E-5 0311 Marine Corps enlisted population group time series introduced in Chapter 
IV, the forecasting model is developed over a multiple-staged process. The first section 
describes the model’s process of smoothing random variation with moving average 
calculations. The following sections demonstrate the procedures for computing 
seasonality and estimating the trend. The chapter concludes by developing the forecast 
model, and creating the E-5 0311 manpower inventory forecast. 

B. MOVING AVERAGES 

Time series decomposition first requires the analyst to smooth the irregular 
component by calculating the series’s moving averages. The forecaster is afforded some 
flexibility in choice of moving average techniques. For time series with higher levels of 
irregularity, Yaffee (2000) recommends longer period moving average calculations. 


Quarterly data should be smoothed with a five-period moving average (Shiskin, Young, 


2 Decomposition can be either a multiplicative or additive model. The multiplicative model is more 
commonly applied (Yaffee, 2000). This thesis describes the multiplicative forecasting method. In the 
additive model, multiplication is replaced by addition and division by subtraction. 
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& Musgrave, 1965). By convention, these are centered moving averages and may be 
weighted or unweighted depending on the analyst’s preference (Stanford, 2002). The 
problem of loss of data at the end of a centered moving average is mitigated by allowing 
the number of elements in the set to diminish to the number of observations remaining as 
the end of the series approaches. This study developed and evaluated numerous 
combinations of moving average period length and weight. The five-month, centered 
moving average described in Equation (6.2) provides the best fit for these particular 
manpower inventory forecasting problems. Equations (6.3) and (6.4) describe how this 
study’s model handles the centered moving average’s loss of data for the final two 
observations in the series. Table 7 provides the moving average calculations for the E-5 


0311 forecasting problem. 











MA, = Dies + Xj By x; + Xi a Ny) (6.2) 
3 

MA,, = 2 %,_3 +%_p +H +H) (6.3) 
4 

MA = D(X,» + X41 +%,) (6.4) 


3 


Period [x [MA__[Period [x [MA__[Period [x [MA 
P1835] 
| T1840] 


| 6] 1857] 1858.4 
| 8] 1877] 1862.2! 
| 1866] 1860.6) 





Table 7. Moving Average Calculations; {,,x5,..., X36} 


This process smoothes the series of the irregular and seasonal components. Thus, 
the values produced by the moving average can be considered a product of each 
observation’s trend and cyclical component (Shiskin et al., 1965). Therefore, Equation 


(6.1) is modified in the following manner: 
a2 


MA=T*C 
x=MA*S*1 (6.5) 


.—=S$*] 
MA 


The ratio x/MA is referred to as the actual-to-moving-average ratio. This 
calculation allows the analyst to isolate the series’s seasonal and irregular components. 
Table 8 depicts the actual-to-moving average ratio calculations for the forecasting 


problem example. 


GEL 
a oS | 13] 0.945266] 25] 1.003495 
-_4| ___ S| 26] 1.001077, 
[____3[ 4.000709 #5] 4.077195[ _27| 0.980041 


6 0.999247} ta] 1.013346[ 30] 1.07363 
[_8[ 1.007948] 20] 1.077839 32 1.014706 
[3 1.002902] 2] 1.045918| 33] 0.993789 





Table 8. Actual-to-Moving Average (x/MA) Calculations; {x,,%,,...,%3¢} 


C; SEASONALITY 

The analyst should now isolate the series’s seasonality component, S. This is 
accomplished via a two-step process. As depicted in Table 9 below, the estimated 
seasonal index for each month is calculated by first averaging all the ratios for that 
particular month across all years of the observation set, and then, if required, 
renormalizing the ratios so that they sum to exactly the number of periods in a season. 
As the total of the monthly averages in Table 9 do not sum to twelve, the indices are 
adjusted in this particular example by multiplying each by the ratio of 
12.00000/11.98964. The sum of the indices should correspond to the number of periods 
in a year; e.g., 52 for time series with weekly observations, twelve for monthly 


observations, and four for quarterly observations. The right-hand column of Table 9 
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depicts the time series’s isolated seasonal component, S$, which represents the percentage 


of normal typically observed in a particular month or season. 


Adjusted 
FY 2002 |FY 2003 |FY 2003 Index 


Oct | (0.945 266] 1.003495] 0.97438] 0.975222 
Nov | i: 1.011665] 1.001077] 1.006371] 1.007241 
Dec___| 1.000109] 1.017195] 0.980041] 0.999115] 0.999978 
Feb | 0.991025] 1.010426] 1.012247] 1.004566] 1.005434 
Mar _|_0.999247| 1.013346] 1.01363] 1.008741] 1.009612 


Apr 1.018747| 1.001766] 0.997976| 1.006163] 1.007032] 
fAug | 1.011599] 0.967117] 1.009724] 0.996147| 0.997007] 
Sept__| 1.007516] 0.995778] 0.993331 0.998875] 0.999736] 


Sum 11.98964 12 





Table 9. Time Series Decomposition E-5 0311 Seasonal Index Calculations 


The forecaster can now de-season the data. This is accomplished by dividing 
each value of the time series by the seasonal index, §. Thus, dividing Equation (6.1) by $ 


yields: 
aia cae (6.6) 


The seasonally adjusted data, x/S, is depicted in the right-hand column of Table 
10 below. In this example, the adjusted seasonal indices, S, are all very close in value to 
1.0 suggesting the seasonality component has little effect on this series. 
D. TREND 

After calculating and isolating the time series’s seasonality, the analyst can now 
estimate the trend component, 7. This is accomplished by regressing the de-seasoned 
inventory level on the period number. Equation (6.7) provides the linear regression 
model calculated for the E-5 0311 forecasting problem. Table 11 depicts the trend 


calculation, 7, for each observation. 


T =1895.075 + 0.7873* period (6.7) 
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Moving Adjusted |Deseasoned 
Inventory |Average ee Data 
Observation |x MA x/MA x/S 


es 
| 1840 i Progra 
-—~3|__1833| 1832.8] 7.000709] 0.999978| 1833 
-——~s|__1857| 1858.4] 0.999247| 1.009612| 1839 
- a err] rasea| oorese| 1 017s] 
sa at —ta80 2] srr7a2[-0snrts| i825 
[___1i]___ta9[__1@27.8| 1.071599] 0.997007| 1855 


Table 10. De-seasoned E-5 0311 Inventory Calculations 
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Moving Adjusted |Deseasoned |Least Sq 
Inventory | Average an Data soar 
Observation |x MA x/MA ss 


eer 
2 1840} r Poor 
/ 31833] 1832.8] 1.000709] 0.999978] 1833] 1897.4368| 
-——s|-_1857| 1858.4 0.999247] 1.009612] 1839] 1899. 7985| 
al err rasa oorede| 1 ovata] 1st 1901 3731 
sa} ast e803 0.8780 0 saorts| Teas} T0276 
[___11|_ 1849] 1827.8] 1.071599] 0.997007] 1855] 1908.7348| 





Table 11. Time Series Decomposition E-5 0311 Trend Calculations 


E. MULTIPLICATIVE DECOMPOSITION FORECASTING 

After computing the adjusted seasonal indices and trend calculations, the analyst 
can construct a forecasting model. Equation (6.8) provides the general multiplicative 
decomposition-forecasting model, while Equation (6.9) details the specific model 


applicable to the E-5 0311 forecasting problem. 
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& =S.*T (6.8) 


x, = S, *(1895.075 +.0.7873i) (6.9) 


Table 12 presents the predicted inventory level for each forecast period, the 
observed (actual) inventory level, and the absolute percentage error. The absolute 
percentage error, described in Equation (5.25), is one of many possible measures of 
forecast effectiveness. Figure 7 provides a graphical presentation of the observed and 


forecast inventory levels for each period. 


Multiplicative Decomposition Forecast Results 
Forecast Period |Actual Inventory|Predicted Inventory |Absolute Percentage Error 


6-Month 1869 1947 4.17% 
12-Month 1880 1932 2.77% 


Table 12. Multiplicative Decomposition Forecast Results 
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Figure 6. Multiplicative Decomposition Forecast Results: E-5 0311: Oct 01-Sept 05 
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F. CHAPTER SUMMARY 

Time series decomposition forecasting is based on the concept that underlying 
factors of a time series can be identified and isolated. These components are trend, cycle, 
seasonality, and irregular occurrences. The analyst develops a forecast model by first 
identifying and removing the component effects from the data (decomposition). After the 
effects are identified and isolated, the analyst creates a prediction by reassembling the 
components (recomposition). This technique is a common and effective method of 


forecasting time series behavior. 
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VII. BOX-JENKINS FORECASTING 


A. OVERVIEW 

In 1970, George Box and Gwilym Jenkins devised a methodology for analyzing 
and forecasting univariate time series data (Hoff, 1983). Unlike the automatic forecasting 
approaches detailed in Chapters V and VI, the Box-Jenkins method is not a forecasting 
algorithm per se. It is a “strategy” for identifying and selecting a forecasting technique 
from a general class of statistical models (Chatfield, 1996). This approach creates 
forecasts by combining two common time series modeling methods, the autoregression 
(AR) and moving average (MA) techniques. The Box-Jenkins process assumes that any 
given time series observation can be modeled as a function of its past values (derived 
through autoregression) and the current and past values of random errors or white noise 


(derived through a moving average). 


This chapter’s purpose is to introduce the concepts of selecting and developing a 
forecasting model with the Box-Jenkins techniques.3 The first section discusses the data 
preparations required to ready the time series for analysis and modeling with the Box- 
Jenkins approach. The following section examines procedures for autoregressive and 
moving average model selection. Next, the chapter reviews the model evaluation 
process. The chapter concludes by discussing the Box-Jenkins forecast created for the E- 
5 0311 Marine Corps enlisted population group introduced in Chapter IV. 

B. DATA PREPARATION 

The Box-Jenkins procedure makes three fundamental assumptions regarding the 
data’s characteristics. First, the procedure requires that the time series be discrete and the 
observations equally spaced over time (Yaffee, 2000). Second, the data must have no 
missing observations. Yaffee (2000) discusses several options for developing algorithms 
to replace missing values if required. Finally, the data must be rendered stationary in the 


mean and variance (Chatfield, 1996). 


3 The scope of this thesis limits the discussion of the Box-Jenkins approach to time series analysis and 
forecasting to a brief overview. For a detailed description of the Box-Jenkins approach to time series 
analysis and forecasting, see Yaffee (2000) and Gujarati (2003). 
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Data preparation entails transformations and differencing. Transformations, such 
as logarithms and square roots, can stabilize variance in an observation set where the 
variation changes with the level. Differencing involves manipulating the time series until 
there are no discernable indications of trend or seasonality. This procedure requires 
taking the difference between sequential observations. Most time series are nonstationary 
and need some differencing and/or transformation prior to Box-Jenkins analysis and 
modeling. The E-5 0311 monthly inventory level time series requires logarithmic 
transformation and first order differencing in order to achieve stationarity. 

1. Transformation— Stabilizing the Variance 

A stationary time series has a constant underlying variance. Often, time series 
data will demonstrate fluctuating variance. Graphing the observations will frequently 
expose this volatility. According to Yaffee (2000), “If the variation in the series 
expands, contracts, or fluctuates with the passage of time, the change in variation will 


usually be apparent in a time plot.” 


If instability in variance is present, the forecaster can affect stationarity via a 
number of possible transformations. A power, logarithmic or Box-Cox transformation 
may stabilize the variance. Power transformations include the square, square root, cube, 


and cube root of the original series. 


A logarithmic transformation effectively stabilizes the variance in the mean of the 
E-5 0311 population group observation set. Table 13 lists the log transformed time series 
observations. 

Z Differencing — Stabilizing the Mean 

Box-Jenkins analysis and modeling requires not only variance stationarity, but 
mean stationarity as well (Yaffee, 2000). This implies no trend or seasonality. 
Deviations about the mean are temporary and the observations demonstrate equilibrium 
about the mean over the term of the series (Sherry, 1984). The analyst should conduct a 
preliminary analysis of the data by simply sketching a time plot and inspecting it for any 
observable patterns in the data. Any clear upward or downward trend or seasonality 
would indicate that the data should be differenced prior to continuing the Box-Jenkins 


analysis. 
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First order differencing is accomplished with the equation: 


Z,=X,-X 


i i i-1? 


ee (7.1) 


where Z; replaces x; in the differenced series. The difference process should be repeated 


as required to achieve stationarity. 


Correlograms are also a useful diagnostic tool for determining the need for 
differencing. Rapid attenuation of the plotted autocorrelation function suggests the series 
is sufficiently differenced (Chatfield, 1996). If the autocorrelations increase, decay 
slowly, exhibit a wave-like cyclical pattern, or decrease linearly, passing through zero to 
become negative, the series is not stationary and should be differenced one or more times 
(Yaffee, 2000). Figure 8 depicts a correlogram of the E-5 0311 time series introduced in 
Chapter IV. The decreasing linear pattern depicted by the shaded graph indicates the 
series is not stationary. Figure 9 depicts the series after first order differencing and 
reflects a series that is sufficiently stationary in the mean to continue the Box-Jenkins 
analysis. The analyst can employ the Dickey-Fuller test if a more sophisticated test for 
stationarity is required (Yaffee, 2000). Table 13 provides the log transformed, first 


differenced E-5 0311 inventory time series observations. 
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Table 13. The log transformed, first differenced E-5 0311 inventory time series 
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E-5 0311 Autocorrelation Functions: No Differencing 


























Figure 7. E-5 0311 Autocorrelation Functions: No Differencing 





E-5 0311 Autocorrelation Function: First Differencing 
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Figure 8. E-5 0311 Autocorrelation Functions: First Differencing 


C; MODEL SELECTION 

The Box-Jenkins approach provides a methodology for selecting, developing, and 
evaluating forecasting models from a broad class of time series processes. These include 
autoregressive (AR) and moving average (MA) techniques. AR models assume that the 


current value of a series is dependent on previous observations with some random error 
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component (Chatfield, 1996). MA models assume the current value of the series is 
dependent on external influences. A combination of these techniques is called an 
autoregressive moving average model (ARMA). This method isolates and models both 
the time series’s underlying process and external influences. This systematic approach 
and the large class of available techniques empower the Box-Jenkins procedure to 
effectively model a wide spectrum of time series behavior. 

1. Autoregressive Models 

The analyst should consider the AR process when it is plausible to assume that the 
current observation of the series is dependent on recent previous values with some 
additional random error (Chatfield, 1996). The Marine Corps enlisted inventory time 
series discussed in this thesis are examples of data that fit this description. With the 
manpower inventory process, the current population level is dependent on the previous 


level. The error term accounts for promotions, discharges, and other losses and gains. 


The error term is subject to several assumptions. It must have a Normal 
distribution, a mean of zero, and a constant variance in order for AR modeling to be 


indicated (Yaffee, 2000). An autoregressive process of order p [AR(p)], is of the form 
Z, = 92Z,,+9Z,5++.4+9,Z,_,+&, Ci2) 
where ¢g denotes the autoregressive coefficients and ¢ is the error term. An AR (1) 
model is defined by Equation (7.3) 
Z,=9Z2,,+ 6, (7.3) 


If the model does not satisfy Equation (7.4) below, the series does not meet 
stationarity requirement and requires additional differencing or transformation (Yaffee, 


2000). 


<1 (7.4) 
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The autocorrelation function (ACF) and partial autocorrelation function (PACF) 


are helpful tools for AR model selection.4. The construction of correlograms simplifies 


the interpretation of ACF and PACF parameters. Table 14 details the AR model 


selection criteria given an ACF and PACF correlogram for a given time series. Figures 


10 and 11 illustrate the autocorrelation and partial autocorrelation behavior of AR (1) and 


AR (2) models, respectively. 


AR Model Selection Criteria 
1. If none of the simple autocorrelations is significantly different from zero, 
the series is essentially a random number or white-noise series, which is not 
amenable to autoregressive modeling. 


oe If the simple autocorrelations decrease linearly, passing through zero to 
become negative, or if the simple autocorrelations exhibit a wave-like cyclical 
pattern, passing through zero several times, the series is not stationary; it must be 
differenced one or more times before it may be modeled with an autoregressive 
process. 


3. If the simple autocorrelations exhibit seasonality; i.e., there are 
autocorrelation peaks every dozen or so (in monthly data) lags, the series is not 
stationary; it must be differenced with a gap approximately equal to the seasonal 
interval before further modeling. 


4, If the simple autocorrelations decrease exponentially but approach zero 
gradually, while the partial autocorrelations are significantly non-zero through 
some small number of lags beyond which they are not significantly different 
from zero, the series should be modeled with an autoregressive process. 


a: If the partial autocorrelations decrease exponentially but approach zero 
gradually, while the simple autocorrelations are significantly non-zero through 
some small number of lags beyond which they are not significantly different 
from zero, the series should be modeled with a moving average process. 


6. If the partial and simple autocorrelations both converge upon zero for 
successively longer lags, but neither actually reaches zero after any particular 
lag, the series may be modeled by a combination of autoregressive and moving 
average process. 





Table 14. Autoregressive Model Selection Criteria (From: Arsham, 1994) 


4 “The partial autocorrelation measures correlation between time series observations that are k time 


periods apart after controlling for correlations at intermediate lags (i.e. lags less than &). 


Partial 


autocorrelation is the correlation between Z, and Z,., after removing the effects of the intermediate Z’s.” 


(Gujarati, 2003) 
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Figure 9. | AR (1) Model ACF and PACF Behavior (From: Borhers & Wessa, 2000) 
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Figure 10. AR (2) Model ACF and PACF Behavior (From: Borhers & Wessa, 2000) 
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2. Moving Average Models 

A time series is influenced by a moving average process if the current observation 
can be expressed as a linear function of the current error and one or more previous error 
terms (Sherry, 1984). This process is one affected by a variety of “random,” external 
influences. These events have an immediate effect and may also affect observation 
values, to a certain extent, for several subsequent periods (Chatfield, 1996). Examples 


would be the effects of strikes or natural disasters on a time series of economic indicators. 


A moving average process of order g [MA(q)]is expressed: 





Ce 4 Oe Goan Oe (7.5) 


q t-q 


where @, denotes the moving average parameters and ¢, is the random error term. This 


error term is subject to the same assumptions as the AR error term. Specifically, it must 
be normally distributed, have a mean of zero, and have a constant variance. Equation 


(7.6) expresses a MA (1) model. 
Z, =é,-66,, (7.6) 


As with the AR model identification and selection process, correlograms plotting 
the ACF and PACF are useful for MA model selection. If the ACF of the differenced 
series displays an abrupt drop to zero while the PACF decays more slowly, the series is 
likely influenced by a MA process. Examination of the ACF plot reveals the likely order, 
q, of the moving average model. If the autocorrelation is significant at lag k, but not at 
any higher lags (the ACF drops to zero at lag k+/), a moving average model of order k 
[MA (k)] is indicated. Figure 12 graphically depicts the behavior of the ACF and PACF 
in MA (1) models while Figure 13 depicts the behavior of those functions in MA (2) 


models. 
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Figure 11. . MA (1) Model ACF and PACF Behavior (From: Borhers & Wessa, 2000) 




















Figure 12. MA (2) Model ACF and PACF Behavior (From: Borhers & Wessa, 2000) 
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Js Autoregressive Moving Average Models 

Autoregressive moving average (ARMA) models integrate the AR and MA 
processes. These models express the relationship of the current observation as a linear 
function of past observations and error terms. As previously noted, an autoregressive 
process of order p is typically denoted as AR (p). A moving average process with q lag 
terms is conventionally classified as MA (q). A combination model containing p AR 


terms and g moving average terms is classified as ARMA (p,q). 


If the time series is differenced d times to meet the constraints of stationarity, the 
model is classified as a ARIMA (p,d,q) model. The “I’ denotes that the model is 


“integrated.” The equation for this model is: 


Z,=O4 OZ, +-4+9,Z,_, +E, -O6,.-- OE (7.7) 


q t-q 


The autocorrelation and partial autocorrelation functions are used to select model 
candidates as described in the previous sections. 
D. MODEL EVALUATION 

The next step in the Box-Jenkins process is to estimate the parameters and 
evaluate one or more models. Parameter estimation involves calculating the model 
coefficients that provide the best fit for the data and is frequently accomplished with 
statistical software. Once the models’ parameters have been estimated, the analyst must 
evaluate them for accuracy assessing the autocorrelation function of the residuals, the t- 
ratios, and computing the minimum sum of squares or some other goodness-of-fit 


indicator for comparison purposes (Sherry, 1984). 


The autocorrelation function of each model’s residuals should be random about 
the mean and the mean should be zero. Additionally, the ACF should have a constant 
variance and a magnitude of less than two standard errors. The t-ratio is calculated for 
each parameter estimate by the series’s standard deviation. This ratio should be greater 
than +/- 2.0 indicating that the coefficient is significantly different from zero. The sum of 
squares measurement is useful for comparing the goodness-of-fit between two or more 


models. 
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Often more than one model will satisfy the evaluation criteria. Under these 
circumstances, the analyst should apply the principle of parsimony and select the lowest- 
order model available that satisfies the analyst’s forecast criteria.5 
E. FORECASTING 

A Box-Jenkins analysis of the E-5 0311 time series data indicates that a log- 
transformed ARIMA (1,1,1) model is appropriate for this particular population group. 
Modifying Equation (7.7), an ARIMA (1,1,1) model can be expressed mathematically as: 


Z,= A, + QZ, + &,— AE, (7.8) 


Recalling Equation (7.1), this formula is rewritten to the following form for 
models with first order differencing (O’Donovon, 1983): 
X,— X41 =O +O (% 4 -%2) +, - GE 


(7.9) 
X= A, dice J), = AX,» + €,— HE,., 


Table 15 details the model’s parameter estimates calculated by the SAS Time 


Series Forecasting System software application. 


Model Parameter [Estimate |StdError|t = —‘[Prob>|t| 
-0.0003106] 0.0058] -0.0534] 0.9577, 


Model Variance | o.ooo7ezal [| st ———sY 


Table 15. E-5 0311 Log ARIMA (1,1,1) Parameter Estimates 
Substituting the ARIMA process coefficients outlined in Table 15 above, the 





analyst can mathematically express the E-5 0311 inventory data as:6 
Z, = —0.0003 106 —0.5296*Z,_, +¢,+0.94190e, _, (7.10) 


Table 16 presents the predicted inventory level for each forecast period, the 
observed (actual) inventory level, and the absolute percentage error. The absolute 
percentage error, described in Equation (5.25), is one of many possible measures of 
forecast effectiveness. Figure 14 provides a graphical presentation of the observed and 


forecast inventory levels for each period. 


5 Chapter II discusses forecast criteria. 


6 This model requires a log to linear conversion to complete the inventory forecast computations. 
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Log ARIMA (1,1,1) Forecast Results 


Forecast Period |Actual Inventory|Predicted Inventory |Absolute Percentage Error 


1869 1830 2.09% 
12-Month 1880 1833 1.44% 


Table 16. E-5 0311 Log ARIMA (1,1,1) Forecast Results 
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Figure 13. Log ARIMA (1,1,1) Forecast Results: E-5 0311 -— Oct 01-Sept 05 


The Box-Jenkins technique is not an automated forecasting method. Each time 
series model requires separate analysis, development and evaluation. The log ARIMA 
(1,1,1) model developed for the E-5 0311 time series is not suitable for every Marine 
Corps enlisted population group considered in this thesis. The Appendix details the Box- 


Jenkins forecast developed for each series. 


Because the technique emphasizes the analysis and evaluation processes, it 
frequently yields more effective models than generically applying automated forecasting 
methods. However, the Box-Jenkins method is a complex approach that often requires a 
“large expenditure of time and effort” (Chatfield, 1996). Although statistical computer 


applications can facilitate the development of Box-Jenkins models, automatic techniques 
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are generally more suitable when there are “large numbers of series to be forecast” 
(Chatfield, 1996). 
F. CHAPTER SUMMARY 

The Box-Jenkins technique is a sophisticated approach by which to analyze time 
series data and extrapolate a forecast. This method provides a framework for identifying 
and selecting a model and creates forecasts by combining two common time series 
modeling methods, the autoregression and moving average techniques. These processes 
assume any given time series observation can be modeled as a function of its past values 
and the current and past values of random errors. The Box-Jenkins technique is 
appropriate for developing forecasts for a small number of series or when the analyst has 
a significant interest in diagnosing and analyzing the fundamental components affecting a 


particular time series. 
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VUI. RESULTS 


A. OVERVIEW 

This chapter evaluates the six- and twelve-month predictive capability of the time 
series forecasting models developed in this thesis. The first section provides an overview 
of evaluation techniques and considerations. Several useful measures of effectiveness are 
introduced. The following sections compare the forecast results with respect to the Mean 
Absolute Percentage Error, the Sum of Squared Errors, and the Mean Absolute Error 
criteria respectively. The chapter concludes with an analysis of the results. 
B. EVALUATING FORECASTS 

Forecast evaluation is an essential part of the model development process. 
Methods for evaluating forecast models are based on conventional analytical procedures. 
The analyst should evaluate the models under situations that match the forecasting 
problem. One effective method is to back-forecast, employing older observations in the 
time series to develop the model, while comparing the forecast fit against the actual 
values of recent, out-of-sample observations (Wooldridge, 2003). Yaffee (2000), notes 
that a model’s statistical accuracy should not be the single criterion employed when 
evaluating a particular forecasting technique. The analyst should also develop and 
consider some assessment of cost in terms of time, money, and effort required for each 


method. 


When evaluating the effectiveness of forecast methods across several series, the 
analyst should utilize error measurements that consider the scale of the observations 
(Armstrong, 2001). This measurement must not be subject to distortion due to a disparity 
in series observation value magnitudes (thus weighting some samples more heavily). 
Some forecasts perform well when evaluated against one criterion yet are outperformed 
by other models when judged by other criteria (Yaffee, 2000). Consequently, the analyst 
should employ and consider multiple error measures when evaluating models 


(Armstrong, 2001). 


This thesis develops manpower inventory forecasting models for 44 


MOS/paygrade combinations. The Appendix summarizes the time series forecasts 
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constructed for each population group. Each observation set contains 48 monthly 
observations starting in October 2001 and concluding in September 2005. The estimation 
sub-sample consists of the first 36 observations and is used to estimate each model’s 
parameters. The validation sub-sample consists of the final 12 observations and is 
reserved for evaluating each model’s forecast accuracy with respect to actual inventory 


levels. 


Three error measurements are offered to comparatively evaluate the models 
against each other and the current forecasting method.?7 These are the Mean Absolute 
Percentage Error (MAPE), Sum of Squared Errors (SSE), and Mean Absolute Error 
(MAE). Equations (8.1), (8.2), and (8.3) respectively describe these assessment-of-fit 


indicators. 





MAPE = #100 (8.1) 

T 
SSE =) (%,-%&Y (8.2) 

i=l 

T 

DIG ~X;) 
MAE = =1__ (8.3) 
c 


where 7 = number of forecasts analyzed, 
i = the forecast period, 
x, = the actual observed value at time i, 


and x, =the forecast value at time i. 


C. RESULTS 

1. Mean Absolute Percentage Error 

The Mean Absolute Percentage Error is a general assessment of fit useful in 
comparing the effectiveness of different models (Yaffee, 2001). The scale of the time 
series observations influences this particular error measurement. Consequently, a 


forecast error of 1,000 in a population group of 10,000 Marines carries the identical 


7 The Enlisted Staffing Goal Model’s inventory forecasting method, described in Chapter III, produces 
a single forecast — regardless of the forecast horizon. Consequently, in this evaluation, the current 
method’s six- and twelve-month forecasts are the same. 
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weight as a forecast error of 10 in a population group of 100 Marines. Table 17 provides 


the six- and twelve-month MAPE comparisons for the forecasting techniques considered 


in this study. A lower value indicates a better forecast fit. 





Six-Month Forecast 


Twelve-Month Forecast 























Current Method 10.2% 10.7% 
Exponential Smoothing 11.0% 12.7% 
Decomposition 12.3% 13.2% 
Box-Jenkins 10.5% 10.7% 








Table 17. 


2 Sum of Squared Errors 


Forecast Model Comparison - Mean Absolute Percentage Error (MAPE) 


The Sum of Squared Errors is a common statistical method for evaluating and 


comparing a model’s accuracy. 


This statistic provides the analyst a sense of the 


forecast’s dispersion of error (Yaffee, 2000). Table 18 compares the respective models’ 


six- and twelve-month forecast SSE calculations. A lower value indicates a superior fit. 


Current Method 
Exponential Smoothing 
Decomposition 


Box-Jenkins 





Six-Month Forecast 
1,529,234 
475,752 
196,944 
475,045 





Twelve-Month Forecast 
1,236,988 
300,874 
159,051 
248,693 





Table 18. 


3. Mean Absolute Error 


Forecast Model Comparison - Sum of Squared Errors (SSE) 


The Mean Absolute Error is a weighted average of the absolute errors, with the 


relative frequencies as the weight factors. 


Simply put, it represents each forecasting 


technique’s average error amount. Table 19 compares the models’ six- and twelve-month 


forecast SSE calculations. A lower value indicates a superior forecast fit. 
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Six-Month Forecast 


Twelve-Month Forecast 























Current Method 47 49 
Exponential Smoothing 34 37 
Decomposition 28 30 
Box-Jenkins 30 32 





Table 19. 


4. E1-E4 Paygrade Forecast Evaluation 


Forecast Model Comparison — Mean Absolute Error (MAE) 


Table 20 lists the six- and twelve-month forecast error measurement calculations 


for a sub-set of the population group forecasts included in this study. 


comprised of the ten E1-E4 paygrade time series. 


This sub-set is 























6-Month Forecast 12-Month Forecast 
MAPE SSE MAE MAPE SSE MAE 

Current 

Method 13.4% | 1,508,894 38 15.9% 1,214,237 39 
Exponential 

Smoothing 10.3% 455,648 23 17.2% 247,724 22 
Decomposition 9.9% 176,113 16 21.5% 136,510 19 
Box-Jenkins 10.8% 467,770 22 12.3% 225,494 21 























Table 20. E1-E4 Paygrade Time Series Forecast Model Comparisons 


D. ANALYSIS 
On the whole, the accuracy of time series forecasts improves, as more 
observations are included (Yaffee, 2000). This axiom is dependant on the data being 
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relevant to the forecast (Yaffee, 2000). Under normal circumstances many of these 
forecast techniques would benefit from an additional year or two of monthly inventory 
level observations. However, these data sets are constrained to 36-months in length to 
better reflect the impact of policy changes since Operation IRAQI FREEDOM and the 
Global War on Terror. The estimation sub-sample still includes a 6-month period of 
Stop-Loss, an irregular component that had a significant effect on many of the population 
group forecasts. Despite the truncated time series span and the conspicuous effects of the 
irregular component, the study’s forecasting models are generally effective. 
Incorporating larger observation sets after manpower levels stabilize in the future will 


likely yield more accurate forecasts. 


While not the best technique evaluated, the Enlisted Staffing Goal Model’s 
current forecast is surprisingly effective relative to the more sophisticated techniques 
developed in this thesis. The more distant into the future a forecast is made, the less 
accurate it generally is. Notably, the current model provides similar levels of 
effectiveness for both six- and twelve-month forecasts. Therefore, Headquarters Marine 
Corps could extend this method’s prediction horizon from six to twelve months if 


required, with little bearing on the forecast’s accuracy. 


While the current method performs modestly better under the MAPE criterion 
than the other methods, the Holt-Winters exponential smoothing, multiplicative 
decomposition, and Box-Jenkins forecasting techniques outperform the status quo under 
the SSE and MAE goodness-of-fit measurements. The multiplicative decomposition 
method was superior to the other forecasting models using the SSE and MAE evaluation 
criteria. This thesis concludes that the multiplicative decomposition forecasting 
technique is the most effective univariate modeling technique to forecast future Marine 


Corps enlisted personnel inventory levels. 


Of note is the six-month forecast returned by the decomposition method for the 
E1-E4 paygrade time series. For these paygrades and this specific forecast horizon, this 
technique surpasses the others in all evaluation criteria considered. This finding warrants 
further research of the technique’s effectiveness on a larger pool of El-E4 population 


group time series. 


a 


E. CHAPTER SUMMARY 

This chapter evaluates Marine Corps enlisted personnel inventory forecasts 
developed by the Holt-Winters exponential smoothing, multiplicative decomposition, and 
Box-Jenkins techniques against the Enlisted Staffing Goal Model’s current forecasting 
method. The predictions are evaluated against actual, out-of-sample inventory levels 
using several goodness-of-fit indicators including Mean Absolute Error Rate, Mean 
Absolute Error, and Sum of Squared Errors. Among the modeling techniques evaluated, 
the multiplicative decomposition performed the best overall and represents an 


improvement over the Marine Corps’ current naive forecasting method. 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

Accurately forecasting future manpower inventory levels in a fundamental 
prerequisite for the development of an effective and functional staffing plan. The 
Enlisted Staffing Goal Model (ESGM) generates a proposed enlisted-force staffing plan 
that mitigates the dilemma of how to optimally meet service-manning requirements 
within the constraints of the Marine Corps’ limited inventory of assignable personnel. In 
order for the staffing plan to be effective, the model’s inputs must be accurate. One of 
these inputs, the forecast inventory available for assignment, is the weakest component of 
the ESGM process. Developing and incorporating a more accurate manpower inventory 
forecasting model will result in a closer to optimal distribution of personnel to meet the 


Marine Corps’ staffing requirements and improve operational readiness. 


A time series is a collection of observations in time. Although the values of 
individual observations cannot be predicted exactly, the distribution of stochastic time 
series observations commonly follows a discernable pattern. Statistical models can often 
describe these patterns. These models assume that the observations vary randomly about 
an underlying mean value that is a function of time. The time series may also be 
characterized by one or more behavioral components that may be isolated, modeled and 


incorporated into the forecast algorithm. 


Univariate time series forecasting models make predictions by extrapolating the 
past behavior of a single variable of interest. Forecasting is appropriate for stochastic 
time series data when the underlying causes of variation do not change significantly over 
time. The forecasting process consists of five steps. This heuristic method involves 
formulating the problem, obtaining and preparing the data, selecting and applying 
forecasting methods, evaluating models, and using forecasts. Guided by this 


developmental framework, analysts can produce more accurate and efficient forecasts. 


This thesis develops and evaluates manpower inventory forecasting models for 44 
representative Marine Corps enlisted population groups using three univariate time series 


forecasting techniques. These methods are the Holt-Winters exponential smoothing, 
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multiplicative decomposition, and Box-Jenkins autoregressive integrated moving average 
forecasting models. Historical personnel strength levels obtained from the Marine Corps 


Total Force Data Warehouse comprise the observation sets employed in this study. 


The Holt-Winters exponential smoothing method is a common univariate time 
series forecasting technique. This model smoothes irregular fluctuations by using 
weighted averages in an exponentially decreasing manner. The technique creates 
forecasts by estimating and extrapolating a linear trend while adjusting the data in each 


period by estimated seasonal indices. 


Multiplicative decomposition is another effective method of forecasting univariate 
time series data. This technique is based on the concept that the underlying factors can 
be identified and isolated. These factors are trend, cycles, seasonality, and random 
variation. The analyst develops a model by first identifying and removing the component 
effects from the data. After these effects are isolated, the analyst creates a forecast by 


reassembling the modeled components. 


The Box-Jenkins technique is a sophisticated approach by which to analyze time 
series data and extrapolate a forecast. This methodology provides a framework for 
preparing data, selecting a model, and creating forecasts and is commonly referred to as 
autoregressive integrated moving average (ARIMA) forecasting. This technique assumes 
any given time series observation can be modeled as a function of its past values and 


current and past values of random errors. 


Models developed in this study using the Holt-Winters exponential smoothing, 
the multiplicative decomposition, and the Box-Jenkins forecasting techniques are 
compared to the current forecasting method and each other be applying a variety of 
statistical goodness-of-fit measurements. The evaluation techniques used are the Mean 
Absolute Percentage Error, the Sum of Squared Errors, and the Mean Absolute Error. 
Among the forecasting techniques evaluated, the multiplicative decomposition method 
performed best overall and represents an improvement over the Marine Corp’s current 


forecasting method. 
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B. RECOMMENDATIONS 

This thesis recommends Marine Corps Systems Command, Total Force 
Information Technology Systems develop and introduce a multiplicative decomposition 
forecasting model into the Enlisted Staffing Goal Model. This forecasting technique 
should be implemented in phases, starting with the E-1 through E-4 population groups. 
C. SUGGESTED FURTHER STUDIES 

Although there are no rigorous precepts for the number of time series 
observations required for effective forecast model parameter estimation, the literature 
generally suggest that the accuracy of predictions improve as the sample size increases. 
This thesis limits the estimation sub-sample to 36 observations to better reflect the impact 
of policy changes since the commencement of combat operations in Iraq and elsewhere in 
support of Global War on Terror. Incorporating larger sample sizes will likely yield 
more accurate forecasts. This thesis suggests future studies develop and evaluate models 


employing time series with 48 and 60 observations. 


The scope of this study is limited to univariate time series forecasting techniques. 
Other approaches may be applicable to the manpower inventory level-forecasting 
problem. This thesis recommends future studies investigate the merit of techniques such 


as Markov chains and multivariate forecasting approaches. 
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APPENDIX 


This Appendix contains the actual six and twelve month inventory level and 
forecast results for each population group included in this study. Forecasts for the Holt- 
Winters exponential smoothing, multiplicative decomposition, Box-Jenkins, and the 
current method are provided in the tables below. The Absolute Percentage Error 
computation, described in Equation (5.25) is also provided. The tables also denote the 
specific ARIMA model specified for each time series during the Box-Jenkins model 


identification process. 
































Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 130 121 6.9% 
Current Method 

12-Months 146 121 17.1% 
Holt-Winters 6-Months 130 167 28.5% 
Exponential 
Smoothing 12-Months 146 202 38.4% 
Multiplicative 6-Months 130 159 22.3% 
Decompositow 1 o-touths 146 168 15.1% 
Box ienking 6-Months 130 146 12.3% 
Box-Cox ARIMA (1,1,1) | 12-Months 146 159 8.9% 














Table 21. 0211 E-5 Forecast Results 
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Absolute 



















































































Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 117 121 3.4% 
Current Method 

12-Months 134 ia 9.7% 
Holt-Winters 6-Months 117 137 17.1% 
Exponential 
Smoothing 12-Months 134 144 7.5% 

eae 6-Months 117 163 39.3% 
Multiplicative 
Decomposition | 12-Months 134 163 21.6% 
Hox leans 6-Months 117 117 0% 
ARIMA (2.0,0)(1,0,0,) | 12-Months 134 111 17.9% 
Table 22. 0211 E-6 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 86 62 27.9% 
Current Method 

12-Months 83 62 25.3% 
Holt-Winters 6-Months 86 66 23.3% 
Exponential 
Smoothing 12-Months 83 61 26.5% 
Multiplicative 6-Months 86 83 3.5% 
Pecompestion | Months 83 81 2.4% 
Boxclenking 6-Months 86 71 17.4% 
Logistic ARIMA(1,1,1,)_ | 12-Months 83 75 9.6% 

Table 23. 0211 E-7 Forecast Results 
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Absolute 






















































































Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 17 17 0% 
Current Method 

12-Months 16 17 6.3% 
Holt- Winters 6-Months 17 15 11.8% 
Exponential 
Smoothing 12-Months 16 12 25.0% 
Multiplicative 6-Months 17 25 35.3% 
Decomposition: 15 -Menths 16 24 50.0% 
Box lenlans 6-Months 17 16 5.9% 
ARIMA (20,6100) 12-Months 16 1S 6.3% 

Table 24. 0211 E-8 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 9776 8571 12.3% 
Current Method 

12-Months 9662 8571 9.4% 
Holt- Winters 6-Months 9776 9125 6.7% 
Exponential 
Smoothing 12-Months 9662 9205 4.7% 
Multiplicative 6-Months 9776 9393 3.9% 
PecOMpOMmOns lip Months 9662 9872 2.2% 
Boeelenkins 6-Months 9776 9109 6.8% 
Logistic ARIMA (11,1) | 12-Months 9662 9215 4.6% 

Table 25. 0311 E-1 — E-3 Forecast Results 
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Absolute 






















































































Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 3071 2855 7.0% 
Current Method 

12-Months 3140 2855 9.1% 
Holt- Winters 6-Months 3071 3028 1.4% 
Exponential 
Smoothing 12-Months 3140 S129 0.3% 
Multiplicative 6-Months 3071 2915 5.1% 
Decomposition: 715 Months 3140 2847 9.3% 
Rex Jenkins 6-Months 3071 3044 0.9% 
BoE LQr ARMA CLD) 12-Months 3140 3082 1.8% 

Table 26. 0311 E-4 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 1869 1757 6.0% 
Current Method 

12-Months 1880 1757 6.5% 
Holt- Winters 6-Months 1869 1762 5.7% 
Exponential 
Smoothing 12-Months 1880 1667 11.3% 
Multiplicative 6-Months 1869 1947 4.2% 
Pecempositons | o-Menths 1880 1932 2.8% 
BoxJeaking 6-Months 1869 1830 2.1% 
COP Me a) 12-Months 1880 1833 2.5% 

Table 27. 0311 E-5 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type Horizon Inventory Inventory Error 

6-Months 1387 1328 4.3% 
Current Method 

12-Months 1293 1328 2.7% 
Holt-Winters 6-Months 1387 1370 1.2% 
Exponential 
Smoothing 12-Months 1293 1354 4.7% 
Multiplicative 6-Months 1387 1425 2.7% 
Decomposition | 1 Months 1293 1378 6.6% 
Rex Tenlans 6-Months 1387 1394 0.5% 
ARIMA (2,0,0)(1,0,0) 12-Months 1293 1401 7.71% 

Table 28. 0369 E-6 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type Horizon Inventory Inventory Error 

6-Months 711 697 2.0% 
Current Method 

12-Months 680 697 2.5% 
Holt-Winters 6-Months 711 716 0.07 
Exponential 
Smoothing 12-Months 680 122 6.2% 
Multiplicative 6-Months 711 667 6.2% 
Decomposition: | 15: Months 680 669 1.6% 
Rexctesbing 6-Months 711 722 1.5% 
ARIMA (2,0,0)(1,0,0) 12-Months 680 723 5.9% 

Table 29. 0369 E-7 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 178 191 7.3% 
Current Method 

12-Months 180 191 6.1% 
Holt- Winters 6-Months 178 208 16.9% 
Exponential 
Smoothing 12-Months 180 204 13.3% 
Multiplicative 6-Months 178 197 10.7% 
Decomposition: | 15-Menths 180 189 5.0% 
Boe Tealkins 6-Months 178 216 21.3 
ARIMA (20,6100) 12-Months 180 217 20.6% 

Table 30. 0369 E-8 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 64 66 3.1% 
Current Method 

12-Months 71 66 7.0% 
Holt- Winters 6-Months 64 67 4.7% 
Exponential 
Smoothing 12-Months 71 68 4.2% 
Multiplicative 6-Months 64 74 15.6% 
PecOMpOMmOns lip Months 7 75 5.6% 
Bex Jeakine 6-Months 64 ye) 14.1% 
ARIMA:(2,0,0)(1,0,0) 12-Months 71 73 2.8% 

Table 31. 0369 E-9 Forecast Results 
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Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 73 58 20.5% 
Current Method 

12-Months 88 58 34.1% 
Holt-Winters 6-Months Fis) 90 23.3% 
Exponential 
Smoothing 12-Months 88 103 17.0% 
Multiplicative 6-Months 73 76 4.1% 
Decomposition: | 15-Menths 88 70 20.5% 
Box leulans 6-Months a 67 8.2% 
ARIMA2,1,00,15)) 12-Months 88 82 6.8% 

Table 32. 0321 E-1 — E-3 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 274 223 18.6% 
Current Method 

12-Months 305 223 26.9% 
Holt- Winters 6-Months 274 306 11.7% 
Exponential 
Smoothing 12-Months 305 347 13.8% 
Multiplicative 6-Months 274 237 13.5% 
PecOMpOMmOns lip Months 305 247 19.0% 
Bus-Jeakine 6-Months 274 227 17.2% 
Box-Cox (1,1,0,1,1,) | 12-Months 305 225 26.2% 

Table 33. 0321 E-4 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 248 252 1.6% 
Current Method 

12-Months 282 252 10.6% 
Holt- Winters 6-Months 248 256 3.2% 
Exponential 
Smoothing 12-Months 282 258 8.5% 
Multiplicative 6-Months 248 259 4.4% 
Decomposition: 15 -Menths 282 262 7.6% 
Box-Jenkins 6-Months 248 249 0.4% 
Bendee, 12-Months 282 248 12.1% 

Table 34. 0321 E-5 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 141 129 8.5% 
Current Method 

12-Months 133 129 3.0% 
Holt- Winters 6-Months 141 130 7.8% 
Exponential 
Smoothing 12-Months 133 126 5.3% 
Multiplicative 6-Months 141 153 8.5% 
PecOMpOMmOns lip Months 133 151 13.5% 
Baselenking 6-Months 141 132 6.4% 
ARIMA:(2,0,0)(1,0,0) 12-Months 133 130 2.3% 

Table 35. 0321 E-6 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 93 91 2.2% 
Current Method 

12-Months 87 91 4.4% 
Holt- Winters 6-Months 93 91 2.2% 
Exponential 
Smoothing 12-Months 87 90 3.4% 
Multiplicative 6-Months 93 90 3.3% 
Decomposition | 19 ‘Months 87 90 3.4% 
Box Teakins 6-Months 93 97 4.3% 
ARIMA (2,0,0) (1.0.0) | 12-Months 87 98 12.6% 

Table 36. 0321 E-7 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 23 26 13.0% 
Current Method 

12-Months 27 26 3.7% 
Holt- Winters 6-Months 23 25 8.7% 
Exponential 
Smoothing 12-Months 29 oa, 0% 
Multiplicative 6-Months 23 30 30.4% 
PecOMpOMmOns lip Months 27 31 14.8% 
Box-lenking 6-Months 23 26 13.0% 
ARIMA:(2,0,0)(1,0,0) 12-Months Dei 24 11.1% 

Table 37. 0321 E-8 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 10 8 20% 
Current Method 

12-Months 11 8 27.3% 
Holt- Winters 6-Months 10 9 10% 
Exponential 
Smoothing 12-Months 11 9 18.2% 
Multiplicative opr ~ We 
Decomposition: 15 Months ll ll 0% 
Box-Jenkins 6-Months 10 9 10% 
fy || Poutis i 9 18.2% 

Table 38. 0321 E-9 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 83 85 2.4% 
Current Method 

12-Months 114 85 25.4% 
Holt- Winters 6-Months 83 109 31.3% 
Exponential 
Smoothing 12-Months 114 118 3.5% 
Multiplicative 6-Months 83 97 16.7% 
PecOmpommons lip Months 114 104 8.8% 
Box-Jenkins 6-Months 83 105 26.5% 
OED” | 12-Months 114 11 2.6% 

Table 39. 2336 E-5 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type Horizon Inventory Inventory Error 

6-Months 110 100 9.1% 
Current Method 

12-Months 101 100 1.0% 
Holt-Winters 6-Months 110 100 9.1% 
Exponential 
Smoothing 12-Months 101 102 1.0% 
Multiplicative 6-Months 110 100 9.1% 
Decomposition | 1 Months 101 100 1.0% 
Box lenlans 6-Months 110 101 8.2% 
ARIMA (2,0,0)(1,0,0) 12-Months 101 102 1.0% 

Table 40. 2336 E-6 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type Horizon Inventory Inventory Error 

6-Months 71 57 19.7% 
Current Method 

12-Months 56 57 1.8% 
Holt-Winters 6-Months 71 45 36.6% 
Exponential 
Smoothing 12-Months 56 28 50.0% 
Multiplicative 6-Months 71 65 8.5% 
Decomposition: | 15: Months 56 64 14.3% 
Box lenkine 6-Months 71 66 7.0% 
ARIMA (2,0,0)(1,0,0) 12-Months 56 65 16.1% 

Table 41. 2336 E-7 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 22 26 18.2% 
Current Method 

12-Months pie 26 13.0% 
Holt-Winters 6-Months 22 26 18.2% 
Exponential 
Smoothing 12-Months 23 24 4.3% 
Multiplicative 6-Months pap) 30 36.4% 
Decomposition: 15 Months 23 29 26.1% 
Box-Jenkins 6-Months 22 28 27.3% 
Op ARMA | 12-Months 23 28 21.8% 

Table 42. 2336 E-8 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months vi 8 14.3% 
Current Method 

12-Months 9 8 11.1% 
Holt- Winters 6-Months i 8 14.3% 
Exponential 
Smoothing 12-Months 9 8 11.1% 
Multiplicative 6-Months 7 9 28.6% 
Decomposition 12-Months 9 9 0% 
Boxcienkins 6-Months 7 8 14.3% 
ARIMA:(2,0,0)(1,0,0) 12-Months 9 8 11.1% 

Table 43. 2336 E-9 Forecast Results 
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Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 1159 1211 4.5% 
Current Method 

12-Months 1147 120 4.5% 
Holt- Winters 6-Months 1159 1301 12.3% 
Exponential 
Smoothing 12-Months 1147 1308 14.0% 
Multiplicative 6-Months 1159 1110 4.2% 
Decomposition: 15 -Menths 1147 1076 6.2% 
Rex Tenlans 6-Months 1159 1290 11.3% 
ARIMA (20,6100) 12-Months 1147 1331 16.0% 

Table 44. 5811 E-1 — E-3 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 618 556 10.0% 
Current Method 

12-Months 679 556 18.1% 
Holt- Winters 6-Months 618 525 15.0% 
Exponential 
Smoothing 12-Months 679 474 30.2% 
Multiplicative 6-Months 618 630 1.9% 
Pecempositons | o-Menths 679 611 10.0% 
Box-lenking 6-Months 618 573 7.3% 
ARIMA:(2,0,0)(1,0,0) 12-Months 679 555 18.3% 

Table 45. 5811 E-4 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 489 494 1.0% 
Current Method 

12-Months 505 494 2.2% 
Holt- Winters 6-Months 489 478 2.2% 
Exponential 
Smoothing 12-Months 505 505 0% 
Multiplicative 6-Months 489 561 14.7% 
Decomposition: | 15-Menths 505 570 12.9% 
Boe Tealkins 6-Months 489 496 1.4% 
ARIMA (20,6100) 12-Months 505 493 2.4% 

Table 46. 5811 E-5 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 300 303 1.0% 
Current Method 

12-Months 259 303 17.0% 
Holt- Winters 6-Months 300 292 2.7% 
Exponential 
Smoothing 12-Months 259 271 4.6% 
Multiplicative 6-Months 300 334 11.3% 
PecOMpOMmOns lip Months 259 332 28.2% 
Bex Jeakine 6-Months 300 321 7.0% 
ARIMA:(2,0,0)(1,0,0) 12-Months 259 313 20.8% 

Table 47. 5811 E-6 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 136 163 19.9% 
Current Method 

12-Months 154 163 5.8% 
Holt- Winters 6-Months 136 177 30.1% 
Exponential 
Smoothing 12-Months 154 193 25.3% 
Multiplicative 6-Months 136 158 16.2% 
Decomposition: 15 -Menths 154 160 2.6% 
Box-Jenkins 6-Months 136 159 16.9% 
Gon. | PoaMoutiis 154 161 4.5% 

Table 48. 5811 E-7 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 56 46 17.9% 
Current Method 

12-Months 50 46 32.0% 
Holt- Winters 6-Months 56 47 16.1% 
Exponential 
Smoothing 12-Months 50 43 14.0% 
Multiplicative 6-Months 56 57 1.8% 
Pecempositons | o-Menths 50 56 12.0% 
Bex Jeakine 6-Months 56 a 1.8% 
ARIMA:(2,0,0)(1,0,0) 12-Months 50 a7 14.0% 

Table 49. 5811 E-8 Forecast Results 
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Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 20 18 10.0% 
Current Method 

12-Months 18 18 0% 
Holt- Winters 6-Months 20 18 10.0% 
Exponential 
Smoothing 12-Months 18 18 0% 
Multiplicative 6-Months 20 16 20.0% 
Decomposition: 15 -Menths 18 15 16.1% 
Box lenlans 6-Months 20 17 15.0% 
EOS RIM) 12-Months 18 17 5.6% 

Table 50. 5811 E-9 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 49 68 38.8% 
Current Method 

12-Months 86 68 20.9% 
Holt- Winters 6-Months 49 61 24.5% 
Exponential 
Smoothing 12-Months 86 40 S307) 
Multiplicative 6-Months 49 32 34.7% 
PecOmpommons lip Months 86 2 97.7% 
Box-Jenkins 6-Months 49 69 40.8% 
GLb | 12-Months 86 57 33.7% 

Table 51. 6092 E-1 — E-3 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 144 163 13.2% 
Current Method 

12-Months 150 163 8.7% 
Holt- Winters 6-Months 144 139 3.5% 
Exponential 
Smoothing 12-Months 150 124 17.3% 
Multiplicative 6-Months 144 169 17.4% 
Decomposition: 15 -Menths 150 171 14.0% 
Box lenlans 6-Months 144 154 6.9% 
ARIMA (20,6100) 12-Months 150 157 4.7% 

Table 52. 6092 E-4 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 124 115 7.3% 
Current Method 

12-Months 118 115 2.6% 
Holt- Winters 6-Months 124 119 4.0% 
Exponential 
Smoothing 12-Months 118 117 0.8% 
Multiplicative 6-Months 124 126 1.6% 
PecOMpOMmOns lip Months 118 125 5.9% 
Baselenking 6-Months 124 120 3.2% 
ARIMA:(2,0,0)(1,0,0) 12-Months 118 119 0.8% 

Table 53. 6092 E-5 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 61 61 0% 
Current Method 

12-Months 56 61 8.9% 
Holt- Winters 6-Months 61 62 1.6% 
Exponential 
Smoothing 12-Months 56 62 10.7% 
Multiplicative 6-Months 61 62 1.6% 
Decomposition: 15 -Menths 56 60 7.1% 
Box lenlans 6-Months 61 61 0% 
ARIMA (20,6100) 12-Months 56 61 8.9% 

Table 54. 6092 E-6 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 27 25 TA% 
Current Method 

12-Months 28 25 10.7% 
Holt- Winters 6-Months 27 24 11.1% 
Exponential 
Smoothing 12-Months 28 22 21.4% 
Multiplicative 6-Months 27 32 18.5% 
PecOmpommons lip Months 28 32 14.2% 
Box-Jenkins 6-Months 27 a2 18.5% 
Goo. | 2eMionths 28 32 14.2% 

Table 55. 6092 E-7 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 200 162 19.0% 
Current Method 

12-Months 175 162 74% 
Holt- Winters 6-Months 200 161 19.5% 
Exponential 
Smoothing 12-Months 175 159 9.1% 
Multiplicative 6-Months 200 175 12.5% 
Decomposition: 15 -Menths 175 175 0% 
Box-Jenkins 6-Months 200 166 17.0% 
Gon. | PoeMioutlis 175 165 5.7% 

Table 56. 6019 E-8 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 79 81 2.5% 
Current Method 

12-Months 87 81 6.9% 
Holt- Winters 6-Months 79 85 7.6% 
Exponential 
Smoothing 12-Months 87 85 2.3% 
Multiplicative 6-Months 719 87 10.1% 
PecOMpOMmOns lip Months 87 88 1.1% 
Box-lenking 6-Months 79 83 4.8% 
ARIMA:(2,0,0)(1,0,0) 12-Months 87 82 5.170 

Table 57. 6019 E-9 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 213 219 2.8% 
Current Method 

12-Months 189 219 15.9% 
Holt-Winters 6-Months 213 218 2.3% 
Exponential 
Smoothing 12-Months 189 225 19.0% 
Multiplicative 6-Months 213 236 10.8% 
Decomposition: | 15-Menths 189 256 35.4% 
Box lenlans 6-Months 213 203 4.7% 
ARIMA (20,6100) 12-Months 189 191 1.1% 

Table 58. 7257 E-1 — E-3 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 137 128 6.6% 
Current Method 

12-Months 141 128 9.2% 
Holt- Winters 6-Months 137 135 1.5% 
Exponential 
Smoothing 12-Months 141 139 1.4% 
Multiplicative 6-Months 137 139 1.5% 
PecOMpOMmOns lip Months 141 140 0.7% 
Box-lenking 6-Months 137 129 5.8% 
ARIMA:(2,0,0)(1,0,0) 12-Months 141 17 9.9% 

Table 59. 7257 E-4 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 113 102 9.7% 
Current Method 

12-Months 112 102 8.9% 
Holt- Winters 6-Months 113 100 11.5% 
Exponential 
Smoothing 12-Months 112 103 8.0% 
Multiplicative 6-Months 113 106 6.2% 
Decomposition: 15 Months 112 107 4.5% 
Box-Jenkins 6-Months 113 112 0.9% 
Op ARMA | 12-Months 112 112 0% 

Table 60. 7257 E-5 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 86 97 12.8% 
Current Method 

12-Months 81 97 18.6% 
Holt- Winters 6-Months 86 89 3.5% 
Exponential 
Smoothing 12-Months 81 77 4.9% 
Multiplicative 6-Months 86 99 15.1% 
PecOmpommons lip Months 81 97 19.8% 
Box-Jenkins 6-Months 86 106 23.3% 
Goo | 2eMfonthis 81 105 29.6% 

Table 61. 7257 E-6 Forecast Results 
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Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 65 54 16.9% 
Current Method 

12-Months 60 54 10.0% 
Holt-Winters 6-Months 65 65 0% 
Exponential 
Smoothing 12-Months 60 70 14.3% 
Multiplicative 6-Months 65 53 18.5% 
Decomposition: 15 -Menths 60 56 6.7% 
Box-Jenkins 6-Months 65 52 20.0% 
Gon. | PoaMoutiis 60 53 13.2% 

Table 62. 7257 E-7 Forecast Results 
Absolute 

Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 19 23 21.1% 
Current Method 

12-Months 19 23 21.1% 
Holt- Winters 6-Months 19 pA 10.5% 
Exponential 
Smoothing 12-Months 19 20 5.3% 
Multiplicative 6-Months 19 a3 17.4% 
PecOMpOMmOns lip Months 19 23 TAG 
Baselenking 6-Months 19 22 15.8% 
ARIMA:(2,0,0)(1,0,0) 12-Months 19 22 15.8% 

Table 63. 7291 E-8 Forecast Results 


84 











Absolute 



































85 








Forecast Actual Predicted Percentage 
Forecast Type | Horizon Inventory Inventory Error 

6-Months 0% 
Current Method 

12-Months 20.0% 
Holt-Winters 6-Months 0% 
Exponential 
Smoothing 12-Months 20.0% 
Multiplicative o-Mouns Oe 
Decomposiwon: |! 45 Menthis 20.0% 
Box-Jenkins 6-Months 0% 
ARIMA (2,0,0)(1,0,0) 12-Months 20.0% 

Table 64. 7291 E-9 Forecast Results 
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